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Abstract 



In a recent article we have discussed the connections between averages 
of powers of Riemann's (-function on the critical line, and averages of char- 
acteristic polynomials of random matrices. The result for random matrices 
was shown to be universal, i.e. independent of the specific probability dis- 
tribution, and the results were derived for arbitrary moments. This allows 
one to extend the previous results to logarithmic moments, for which we 
derive the explicit universal expressions in random matrix theory. We then 
compare these results to various results and conjectures for ("-functions, and 
the correspondence is again striking. 
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1 Correlation functions of characteristic polyno- 
mials 



We first briefly review the result of a previous paper 0, in which we have investigated 
the average of a product of characteristic polynomials of a random matrix. Let X 
be an M x M random Hermitian matrix. The correlation function of 2K distinct 
characteristic polynomials is defined as 

2K 

F 2K (\i, \2k) = ( II det ( A - - X ))- (!) 

a=l 

The average is taken with the normalized probability distribution 

P(X) = iexp-iVTrV(X), (2) 

where V(X) is a polynomial in X. The simplest case consists of a Gaussian distri- 
bution, and the result is easily worked out by the orthogonal polynomial method, 

P(X) = -^exp-^TrX 2 , (3) 

with 

M = N-K. (4) 

Defining the middle-point 

i 2K 

A=^EA«, (5) 

a=l 

and the density of eigenvalues at this point 

1 



p(X) = ^v/I^P, (6) 



we introduced the scaling variables 



2K 

x a = 27rXp(A)(A a -A),with5> a = 0, (7) 

a=l 



and consider the large N limit, finite x a 's limit. In this limit we have shown [l| that, 
N 2K 

exp -(- £ V(A0)*W(Ai, ■ • • , \2k) = 
1 i=i 



..exp(-NK) X^du a _ A 2 {u ir --,u K ) 



(27riVp(A)) — /ll^ ex P-<L 



1 o=l 



Itlll|=l(«a-Il) 



in which the contours enclose all the xi's. It has also been shown that the result is 
universal in this scaling limit, i.e. independent of the specific polynomial V which 
defines the probability distribution. For K = 1, we have 

exp{-^(V(A 1 ) + V(X 2 ))}F 2 (X 1 , A 2 ) = e -*( 2 7riVp(A)) — (9) 
2 x 

with x = 7riVp(A)(Ai — A2), the familiar sine-kernel. 

When all the Aj's are equal, we obtain the 2i^-th moment of the characteristic 
polynomial : 

exp-(NKV(X))F 2K (X,---,X) = 

(2nNp(X)) f H^-exp-z(^0 . (10) 

This contour integration reduces to a simple determinant and one finds, 

exp -(NKV(X))F 2K (X, • • • , A) = (27rtfp(A))*V™ if (11) 
Let us denote the last factor by 

In our previous paper |l|, we have compared ( JTT| ) with the average of the 2K-th 
moment of the ^-function g, HI, for which it has been conjectured that 



[ T dt\a\ + it)? K ^ lK a K {\ogT) K2 (13) 
Jo 2 



1 r T 

T 



where ax is a number theoretic coefficient given by the product of the prime p, 

p V ^0 ml 
and is the same as in (|T2|). Since we have with ( [TT]) an expression valid for all K's, 
we shall extend this comparison beween the moments of characteristic polynomials 
and that of zeta-functions to non-integer K, as will be explained in the next section. 
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2 Non-integer power moment 



Since the result ((TTJ) is valid for any K we may now consider the analytic continuation 



to non-integer K. This requires to continue the coefficient 7^ in ( |P2~D to non-integer 
values. This will be needed for obtaining the logarithmic moment of the characteristic 
polynomials. The non-integer power moments are also interesting by themselves, 
since there exists equivalent studies of fractional power moments of the Riemann 
(-function [|]. 

The factor 7^ = J[l\/ (K+l) ! may be expressed through an integral representation, 

/•oo P -t fi _ P -Kt\2 

iog , K = -l it - [K >-L—± ] (1B) 

which is easily checked by expanding the integrand in powers of e~*. It is more 
conveniently handled if we take it as 

/•oo p-t fl _ P -Kt\2 

^ = H~L ^- VFf 1 ' (16) 

since it may then be splitted into two parts for a < 1. 

Expanding the integrand in the power of e~*, and integrating over i, we obtain 

00 00 

log 7^ = -T(l-a)[K 2 -J2(n + l) a + 2j2(n + K+l) a -\n + l) 

00 

-Y J (n + 2K + l) a -\n + l)} 

n=0 

= -T(l - a)[K 2 - C(-a) + 2C(-a, + 1) - 2i^C(l - a, AT + 1) 
-C(-a,2ir+l) + 2^C(l-a,2^+l)] (17) 

in which the limit a — > 1 is meant. The generalized zeta-function £(2, a) is given by 
X^o( a + n )~ z , and ((0, a) = | — a ; it has the expansion 

. 2r(l — 2) r . ,7T2;. ^2, cos(2?27ra) 



(27T) 1 - L 2 ^ n 



l-z 



,7T2;, sin(2n7ra). . . 

+ C0S (y)E^iz^]- (is) 

Z n=l ' t 
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It is then easy to obtain various results for non-integer K. For instance the K = 
limit is obtained by expanding in powers of K. The term of order K 2 is 

1 



log 7* ^ -K 2 [r(l-a)-r(3-a)]T 



n=0 yn + iy- 1 
= -K 2 [r(l-a)-r(3-a)C(2-a)] 

= K 2 (l + c) + 0(K 3 ) (19) 

where c is Euler's constant, c = 0.5772... 

It is also interesting to evaluate jk for K — ^, since it is needed for computing 
the first moment of the characteristic polynomial or of the Riemann ^-function. In 
this case, we have 

f°° 1 1 e~* 
log 71 = lim — / (- — 7- tttttt) dt 

1 00 
= lim --r(l - a) - 2 x - a + 2)° - (n + 2) Q - 1 ]r(l - a) 

n=0 

= limf— + 2 1 - Q ((2 1+Q - l)C(-a) - (2 a - 1)C(1 - a)T{\ - a) 

a— >1 4 

= T ^log2 + ^log7r + 3C , (-l) (20) 

since C(0) = -|, C(-l) = -s» and £"=o(-l)> + 2)~ s = 1 + (2 1 "* - l)C(s). This 
leads to 71 is 1.1432.... In the literature on Riemann (^-functions concerning the 
moments fll3]), bounds have been conjectured for < K < 1 ||, §[], and they amount 
for the equivalent of 7x to 

1 2 

< 7K < ^npr-^Tn t7\- (21) 



T(K 2 + 1) - ,A ~ T(K 2 + 2)(2 - A") 

We find that our result of ( |i~5| ) indeed satisfies this bound for < K < 1. For 
instance, the bounds require 1.1033 < 71 < 1.1768, and we have found 71 = 1.1432. 
It is easy to verify this bound by expanding around K = 1 and K = 0, and it does 
support the conjecture (ETJ). 
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3 Moment of the logarithm 

We now consider the case for which all the A a 's are equal, and expand the 2K-th 
moment in powers of K : 



F 2K {X,---,X)e 



-NKV{\)+NK 



([det(X-X)e~^ x ^} 2K ) 

OO fi-\ JS\p AT AT 

E SMPog I det ( A - x )\ ~ -v v W + T ]'>. (22) 

P =o P- z z 



^From (|TlD ) we know the same / as 



oo jy-2p' 

/=£ — [log(27riVp(A)r (23) 

p'=0 P • 

provided we neglect, for the moment, the factor 7^- and set it equal to one. We have 
thus made an analytic continuation from integer K to a real variable K, and we have 
expanded in K. We have seen in the previous section that, for K small, 

lK = l + K 2 (l + c) +0(K 3 ) (24) 

where c is Euler constant. This correction of order K 2 gives only a sub leading term 
in ( f23"D compared to log(27riVp(A), and thus we were justified to neglect it. 
Comparing fl22|) and (p^ ), we find 

AT AT 1 O™! 

([log I det(A - X)\ - -V{X) + -] 2m ) ~ __^[log(27rJVp(A)P (25) 
Since there is the expansion (|23|) contains only even powers of K, we have also, 
([log I det(A - X)\ - jV(\) + ^r +1 > = (26) 



It is instructive to verify this for m = ; taking a derivative of (|26| ) with respect to 
A one should verify that 

m _L_)> = ! ia M 



X-X" J X-a 

^V'(X) (27) 



which is nothing but the saddle point equation of the large N limit which determines 
pW II- The solution of this Riemann-Hilbert problem is expressed through G(z) = 
jj < Tr(^^) >= — P(z)\JQ(z)), in which P and Q are polynomials fixed 

by the requirement that G falls like 1/z at infinity Then Re G(x + ie) = V'(x), Im 
G(x + ie) = —7rp(x) in the large N limit. The r.h.s. of the above equation is indeed 
iVReG(A) = f V'(A). 

Our argument may be applied to the moment of the ^-function as well. The 
average of 2K moment of ( function has been conjectured ||, ||| as 

^ [ T dt\C(\ + it)\ 2K * lK a K (logT) K2 (28) 
1 Jo 2 

Expanding the above equation in powers of K, we obtain along the same lines 
1 r T 1 1 9rr)l 

T I Pog|C(2 +m 2m -^^T ( lo g log TT (29) 
provided ax goes to one in the limit K — > 0. Selberg |7j has derived 

\ / Vlog(C^ + m 2m - ^^(bglogTT (30) 

where ImlogC(x) = argC(x). 

Our result (^) may be compared to (|29|) or ([!(]), in which the density of state 
27riVp(A) is replaced by logT. The characteristic polynomial det(A — X) has zeros on 
the real axis in the complex plane of A. Hence this function corresponds to Riemann's 
£(| + iX) = A(X)Z(X), which has zeros on the real A line, where 

A(\)=n-^°^+n-Y-\)> ( 31 ) 
Z(X) = e ie ((±+iX). (32) 

The analogy between | det(A - X)\ and + iX)\ = |A(A)||C(| + zA)| leads to 
log |C(|+«A)| ~ log | det(A-s)|-log \A(X)\, and it may thus correspond to log | det(A— 
X)|-^ + f. 

From those log-moment results, we find that the distribution function is a normal 
Gaussian distribution. For the ( function, or more generally for the L functions, it 
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is known that log|L(| + it)\ is ditributed like a random variable, when t is large, 
with a Gaussian density |7|, ||, ||. The coefficient of (p5|) , 2m!/(2 2m m!) is equal to 
-^T(m + i), which is identical to the coefficient of the following Gaussian integral, 

oo „ 1 I , 

x 2m e~ ax dx = -^T(m + -)a~ m ~^. (33) 

-oo V 7T 2 



Hence we find that our moment (|25f ) does follow a normal distribution. 

We have considered up to now the case of all the Aj's equal ; in the following we 
shall consider two different A«'s, and it will also appear that a normal distribution 
holds for the logarithmic moments. 

4 Moments at two different points 

The formula for F 2 k in (f|) provides also the correlation for two different values of 
the "energies", Ai and A 2 , if we set A a = Ai = • • • = A^ and A;, = Xr+i = • • • = X 2 k- 
For instance, when K = 2, (l\ — 2, l 2 — 2, and h + l 2 = 2K), we have x\ = x 2 = 

—X3 = —X4 = x. Then, the contour integral in (H) becomes 

fdu^^ e" i( " 1+M2) (wi - u 2 ) 2 _ 1 sin 2 x^ 

J (27r) 2 (ui — x) 2 (ui + x) 2 {u 2 — x) 2 (u 2 + x) 2 2x 2 X 2 

where x = 7riVp(A)(Ai — A 2 ), with A = ~(Ai + \ 2 ). Note that the contour integral 



formula of (|8|) has been derived in Dyson's short distance limit [|1(], |TTJ. We consider, 
within this Dyson limit, the large x case : x >> 1. The leading term forlarge x, 
is easily obtained from (|8]). The subleading terms have an oscillatory behavior, but 
we limit ourselves for simplicity to the leading term . Then the leading behavior 
becomes 

([det(A 1 -X)] 2 ' 1 [det(A 2 -X)]^)e- Wiy ( Al ) +/lJV - jV ' 2y ( A2 ) + ' 2jV ~ * - (2-kN pf^ 2 

jr. 2 \ll H~t2 ) 

(35) 

Expanding in powers of l± and l 2 , we obtain from the l.h.s. of fl35|), 

/ = £ (2<l)P 1 {2h)P2 ([log | det(Ax -X)\- ^V(X 1 ) + [log | det(A 2 - X)\ 
Pl ^=o Pi! P2! 2 2 

N N 

~ ^V(\ 2 ) + -r) (36) 



For the r.h.s. of (|35T), we have 



p=0 Z p=0pi,p 2 ,2p=pi+P2 **■ P1-P2- V £>X 

Hence, we obtain when p% + p 2 is even integer, 

([log | det(A 1 - X)\ - ^V(X 1 ) + |r[log | det(A 2 - X)\ - ^V{\ 2 ) + 
- rl^Mlfiog^)], (38 ) 



v 2 y p! ov v/2a; 

where p = {jp\ + P2) /2. When pi + p 2 is an odd integer, the correlation vanishes since 
there is no corresponding term in the r.h.s.. The above result has been derived for 
large x. In this limit, we have obtained, 

([log I det(A! - X) - y V(Ax) - log I det(A 2 -X)\ + jV(X 2 )} 2p ) 
~ ([log I det(Ai - X) - log I det(A 2 - X)\] 2p ) 

- 2(i)^M(lo g [2^]r-2(i)^M(lo g «]r (39) 
2 pi I pi y lx 

where we have expanded the binomial forms, and used fl38l) . The difference V(Ai) — 

V(X 2 ) gives a subleading term, and it has been neglected. 

In the simple case, p = 1 for (|39|), we have a cross-term with two logarithms. By 

taking the derivatives of this cross term, we obtain 

d Xl d X2 (\og I det(Ai - X)\ log I det(A 2 - X)\) = ( Tt x^X Tt \ 2 ~^X ) (40) 
which is two point Green function. We have 

(Tr^— Tr—!—) = N 2 G 2c (z u z 2 ) + N 2 G(z x )G{z 2 ) (41) 
Z\ — X z 2 — X 



where the connected two-point Green function G 2c (z\, z 2 ) has been foun= d in ]T2|. IB . 
There it has been shown that 

N 2 G 2c { Zl ,z 2 ) = -d Zl d Z2 log[l - G( Zl )G(z 2 )] 

u(zi) - u{z 2 ) 



d Zl d Z2 log[- 



z\ - z 2 



' ;< to ^.: 4 .n„ - 1) (42) 



2(z 1 -z 2 ) 2V [(zf-4)(z|-4)]V2 



where u(z) = \[z + V z 2 — 4] . So indeed the result ( p9|) for p = 1 is consistent 
with the previously known results (f42|) in the large N limit. Note that the result 



(HH) has been derived by taking the large N limit first ; hence it is a smoothed 
correlation function, which neglects all the oscillatory terms. (The connected two- 
point correlation function is p 2c (Ai, A 2 ) = — 4^2[G 2c (Ai + ie, \ 2 + ie) + G 2c (Ai — it, A 2 — 
ie) — G 2c (Xi + it, A 2 — ie) — G 2c (A — ie, A 2 + ie)), and it becomes — l/27r 2 A^ 2 (Ai — A 2 ) 2 
for Ai close to A 2 . This result is obtained by smoothing the oscillatory part, while 
the exact result is p 2c (Ai,A 2 ) ~ — sin 2 x/tt 2 N 2 (\i — A 2 ); by taking the large-N limit 
first, the sin 2 x is replaced by 1/2.). 

We noticed that a similar formula exists for the Riemann ^-function, although it 
deals with the imaginary part of the logarithm of the ^-function, from which a study 
of the variance of the number of zeros has been discussed in the literature [14], [15[ . 
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